Cation Transport in Polymer Electrolytes: A Microscopic Approach 
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A microscopic theory for cation diffusion in polymer electrolytes is presented. Based on a thorough 
analysis of molecular dynamics simulations on PEO with LiBF4 the mechanisms of cation dynamics 
are characterised. Cation jumps between polymer chains can be identified as renewal processes. 
This allows us to obtain an explicit expression for the lithium ion diffusion constant D^i by invoking 
polymer specific properties such as the Rouse dynamics. This extends previous phenomenological 
and numerical approaches. In particular, the chain length dependence of Dl, can be predicted and 
compared with experimental data. This dependence can be fully understood without referring to 
entanglement effects. 

PACS numbers: 61.20.Qg,61.25.Hq,66.10.Ed 
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Transport of cations in complex systems is of ma- 
jor relevance in the field of disordered ion conductors. 
Specifically, polymer electrolytes [3, H, S] , using lithium 
salts, have been intensively studied experimentally and 
theoretically due to their technological relevance. Free 
lithium ions (Li + ), uncomplexed by the anions, are the 
desirable charge carriers in electrolytic applications. A 
theoretical description of ionic dynamics in terms of mi- 
croscopic properties is difficult because the dynamics of 
the cations and the polymer segments occur on the same 
time scale [H, In contrast, the dynamics of ions in 
inorganic systems can be characterized by ionic hops be- 
tween permanent sites, supplied by the immobile network 

The phenomenological dynamic bond-percolation 
model (DBP) 0] considers that the long-range ion trans- 
port is enabled by renewal events which lift the block- 
ages in the ionic pathways. Since the dynamics after 
a renewal event is statistically uncorrelated to its past, 
the resulting cationic diffusion constant, Du, is deter- 
mined by a 2 /6r re „ where a 2 and r ren denote the typical 
mean square displacement (MSD) and the time period 
between two consecutive renewal events, respectively. In 
the DBP the renewal process is attributed to a local 
structural relaxation process governed by the polymer 
dynamics. Another fruitful approach to understand the 
mechanisms of cation dynamics is by means of molecular 
dynamics (MD) simulations [1, [1 El, El El • Cationic 
dynamics can be divided into three important mecha- 
nisms: (Ml) motion along a chain ("intra-chain"), (M2) 
motion together with chain segments, using the chain as 
a vehicle ("segmental"), and (M3) jumps between differ- 
ent chains (" inter-chain" ) ; see also [13| . This is sketched 
in Fig[TJ Based on the insight from MD simulations, 
Borodin and Smith [l3|, E3 have recently formulated a 
microscopic transport model. Employing appropriately 
defined Monte-Carlo moves and implicitly using the con- 
cept of renewal process the Li + dynamics has been re- 
produced. Among other things, they have quantified the 
relevance of the variety of cation transport mechanisms 
that contribute to the macroscopic cation diffusivity D^. 




Figure 1: Time scales: n is the time scale for intra-chain ionic 
motion (Ml), T2 is the relaxation time of the polymer chain 
(related to tr: see text for details) (M2), and tz is the waiting 
time of an ion between two inter-chain jumps (M3). 



Our methodology is founded on both the approaches. 
First, we express a 2 in terms of (Ml) and (M2) by ex- 
ploiting the fact that the Li + dynamics is strongly corre- 
lated to the polymer segmental dynamics, which in turn 
can be separated into statistically uncorrelated center- 
of-mass (c.o.m) dynamics of the polymer chain (zeroth 
order Rouse mode) and its internal dynamics (higher or- 
der Rouse modes) [Hj]. This implies, for the Li + ions, 
Du = Dc.o.m + D m (ti,T2,t 3 ). The time scales n, r 2 
and T3 characterize each of the mechanisms (Ml), (M2) 
and (M3), respectively. Extending Ref. [H, we derive 
analytical formulas in terms of Tj with correlations be- 
tween (Ml) and (M2) being taken into account. The 
time scales are extracted from long MD simulations of a 
model polymer electrolyte system. Second, we explicitly 
show that (M3) can be identified as a renewal process so 
that, indeed, a prediction of the long-time behavior, i.e. 
Du, becomes possible. As a central feature we obtain 
the A-dependence (N: number of monomers per chain) 
of the Tj and thus of Du. Two different TV-regimes are 
obtained in agreement with experimental data. 

Atomistic NVT MD simulations are performed for the 
system poly(ethylene) oxide (PEO) and LiBF4 with a 
concentration of EO:Li=20:l (EO: ether oxygen). The 
two body effective polarizable potential employed is de- 
scribed in Ref. [5(. We have simulated this system for 
different chain lengths (N = 24 and N = 48) and at 
different temperatures (400 K < T < 450 K). The re- 
spective densities have been chosen to set the average 
pressure to values of the order of 1 MPa. This letter dis- 
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Figure 2: Average square variation of the average oxygen in- 
dex of one chain to which a Li + is associated during time 
interval t. 

cusses the results for the N — 48 system at T — 450 K 
unless specified otherwise. 

As has been observed before [1, 11 Ell EH a Li+ ion is, 
most of the time, coordinated to a single polymer chain 
through EO atoms, interrupted by infrequent transitions 
between different chains. A Li + ion which is bound to a 
polymer chain is coordinated to a few (as 5) and mostly 
contiguous oxygen atoms. After serially indexing the 
oxygen atoms of a chain in succession we mark the posi- 
tion of the Li + at the chain by the average index, n(t), 
of the enumerated oxygen atoms that belong to its co- 
ordination sphere. To elucidate (Ml), we determine the 
average-square variation (An(t) 2 ) of the average oxygen 
index under the constraint that the Li + ion is attached 
to the same chain during the time interval of length t. 
The result is shown in Fig[2j To a good approxima- 
tion one observes diffusive dynamics (An(t) 2 ) = 2D\t 
where D± is the intra-chain ionic diffusivity. To account 
for the slight deviations from linear behavior we choose 
Di = (An(r2) 2 )/2r 2 (r 2 defined in the next paragraph). 
For later purposes we define 

n = (N - if/^D, (l) 

where t± is a measure of the time it takes for the lithium 
ion to diffuse from one end to the other end of the poly- 
mer chain. Here we obtain n ~ 150 ns (with D\ from 
Fig®. 

To characterize (M2) we first analyze the polymer dy- 
namics. In Figj3(a) we display the MSD go (t) for an 
average oxygen atom (i.e. all oxygen atoms were consid- 
ered for analysis irrespective of the presence or absence of 
Li + near an oxygen atom), characterizing the dynamics 
of the polymer segments. According to our general pro- 
cedure all MSD-functions are computed relative to the 
polymer c.o.m. It exhibits a Rouse-like behavior [l5|, El| 
for short times go(t) °^ t a with a w 0.6, saturating at 
go(t) ~ i? 2 /3 where R 2 is the mean square end-to-end 
distance of the polymer . We have included the theoreti- 
cal Rouse prediction, obtained via numerical summation 

2R 2 ^v-i i_ ex p(-^) 
9(t/r R ) = ^J2 ( 2 ) 

P =i F 

where tr and p denote Rouse time and mode number, 



respectively, and the sum is calculated over the iV — 1 
eigenmodes. It yields a reasonable description of the ob- 
served MSD, using tr — 19 ns; see also [13] • Qualita- 
tively, it is expected that the oxygen atoms which are 
temporarily bound to a Li + ion ought to be somewhat 
slower due to the decrease in the local degrees of free- 
dom. This was checked by calculating g ound (t) for those 
oxygen atoms which, during the whole time interval of 
length t, are bound to one Li + ion. Indeed, g ound (t) is 
also consistent with the Rouse prediction, using a longer 
Rouse time r 2 = 42 ns. Naturally, t-iJtr > 1 reflects the 
immobilization of the polymer segments due to the ions. 
Note that for longer t less oxygen atoms contribute to 
this curve so that the statistics gets worse. 

Switching to the Li + dynamics, we first calculate the 
MSD g K L [ 2 (t) of Li+ ions for which \n{t) - n(0)| < 1; see 
FigS^a). For t > 2 ns this curve is close to go° und (t). 
Thus, we can conclude that the Li + motion strictly fol- 
lows the oxygen dynamics in the absence of (Ml). In 
other words, the cations and the polymer segments ex- 
hibit coupled dynamics. Additionally, we find that for 
shorter times the Li + ions are slower than the correspond- 
ing oxygen atoms. 

In the absence of ion jump events between chains (M3) 
one would have Dm = 0. We have identified the jumps 
from a microscopic analysis of the trajectories. On aver- 
age after T3 = 110 ns a Li + ion jumps between two chains. 
To characterize the effect of these jumps we have first de- 
termined 5i{ 123 (t) which is the MSD of a Li+ ion between 
times [to — t, to + 1] if at time to a jump happens and dur- 
ing the intervals [to — t, to] and [to, to + t] the ion stays 
with the same chain, respectively. The MSD is averaged 
over all jumps (Figj3^b)) . Furthermore we have deter- 
mined the MSD <7£j±(t) during the intervals [t — t,t ] 
or [to,to + t], i.e. just before or after an inter-chain jump, 
respectively. For symmetry reasons one has g^+{t) = 

ffii 1 -^) = 9iA 12 (t)- If' additionally, the inter-chain jump 
at to serves as a renewal process, the statistical indepen- 
dence of cationic dynamics before and after the jump re- 
quires 5 £{ 123 (t) = g%f(t) + g%*(t) = 2g^(t). In case 
of correlations a smaller factor is expected compared to 
2. As exhibited in Fig(3h this relation is indeed found 
with minor deviations (prefactor 2.2 instead of 2.0), thus 
validating the fact that inter-chain transitions can be re- 
garded as renewal processes. 

In the following, we derive an explicit expression Dm = 
Dm(ti,T2,t 3 ). Based on the observed correlations be- 
tween a Li + ion and the polymer dynamics, g^ 12 (t) can 
be described by taking into account the Rouse dynam- 
ics (M2) plus the additional intra-chain diffusion (Ml). 
Formally, one can write 

g^ 2 {t) = m{t)-mf)MX2 (3) 

The average is over the probability density that the ini- 
tial oxygen index to which an ion is linked is i and at a 
time t later is j and the distribution of monomer displace- 
ments as predicted in the Rouse theory. From Ref.(l5l] 
one obtains, using (cos 2, (pir(j — 0.5) /N)) mi = 1/2 and 
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Figure 3: (a) MSD of (i) all oxygen atoms (solid), (ii) oxygen 
atoms which are bound to one Li + during time duration t 
(dash), (iii) Li + which changed the average index of its oxygen 
neighbors by at most 1 (dot). Also shown are the Rouse 
predictions with tr=19 ns (□) and T2=42 ns (a), (b) MSD 
of Li + under different constraints (see text). 
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Assuming Gaussian dynamics for (Ml), i.e. (An 2 ) oc 
t (correction to the deviation as observed in FigJ5] can 
be easily implemented but are not of relevance here), 
one finds (cos(pn(i — j)/N))mx = cx Y>(—p 2 t/ T i)- Eq. 0] 
simplifies to (using Eqj2j) 



(5) 



with 1/t%2 = 1/ri + 1/t2. Thus, the dynamical effects of 
(Ml) and (M2) appear through the resulting relaxation 
rate I/V12. Finally, using the renewal property one can 
write Dm = a %i/^ T 3 explicitly as 
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(^ 12 (r)>M3/6r 3 



(6) 



where & 2 M corresponds to the average MSD between suc- 
cessive inter-chain hopping events due to (Ml) and (M2). 
The average is taken over the distribution of time inter- 
vals between these events. For the numerical analysis 
(see below) we take a simple exponential distribution. 

Approximate analytical expressions can be obtained 
by converting the sum in Eq|4] into an integral from 
to oo 15]. Then one obtains g^ 12 = 2R 2 e n-^ 2 y/t/T 12 
for t <C T\2 and gf^ 12 = R 2 /3 for t T12, respectively. 
Inserting these results into EqJ5] gives 
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if t 3 < T12 (7a) 



if t 3 > T12 (7b) 



The scaling of Dm with r 3 in Eq. fTaj) is consistent with 
the numerical results, obtained in Ref.[13]. Eq. fTaj) holds 
for long chains and takes into consideration the implicit 



N 


Ti [ns] 


r 2 [ns] 


r 3 [ns] 


48 


150 


42 


110 


24 


34 (36) 


10 (10) 


90 (110) 



Table I: The relevant time scales n for iV = 48 and iV = 24 
as well as a check of the scaling relations. 



correlations of (Ml) and (M2). By neglecting these 
correlations, as done in Ref.[13j, the term U 1 /r 12 = 

y/l/ri + 1 /r 2 would instead become yl/ri + \/1/t2- 
This would largely overestimate Dm (for instance, in 
PEO/L1BF4 by 35 %) and the contribution of (Ml) for 
the case r 3 ri2- 

In contrast to the DBP-model we obtain Dm oc l/T ren 
only for short chains. The main conclusion, however, that 
the Li + diffusion has the same temperature dependence 
as the inverse Rouse time and thus as the polymer dy- 
namics, remains valid because all the three time scales t; 
have a similar temperature dependence (data not shown) . 

In Tab. I we have compiled ti, T2, t 3 , obtained from our 
simulations for TV = 24 and N = 48. Of major impor- 
tance are their scaling properties with TV. One expects 
ti,t 2 cx TV 2 and r 3 oc TV . Considering the appropriate 
number of eigenmodes, the predictions for TV = 24, based 
on TV = 48, are given in parenthesis. The agreement is 
convincing. By setting t\ — > 00 one can estimate the 
contribution of (Ml) to Dm in the limit of long chains 
(t 3 -c T12). It is as small as 12% for the present system. 

Naively one might expect that a 2 M oc TV and thus 
D M oc TV for all TV, see e.g. because (M1)-(M3) are 
related to local motions. Using the scaling results for the 
Ti one obtains, however, Dm oc TV only for long chains 
(t 3 <C T12) whereas Dm oc TV for short chains. The reason 
is that for short chains a 2 M is limited by the end-to-end 
distance of the polymer which brings an TV-dependence. 

In the present case, the crossover in Dm, i-e. t 3 w T12, 
occurs for TV sa 100 (see inset of FigH]). It has been specu- 
lated [l8| that the emergence of the entanglement regime 
(TV w 75 [H, [l!|) leads to the crossover with an accom- 
panying change of cation conduction process from poly- 
mer c.o.m dominated motion to percolation type trans- 
port (i.e. (M3)). However, we find, it is a coincidence 
that the entanglement length is similar to the crossover 
length. Thus, the crossover to Dm oc TV is physically 
unrelated to entanglement effects. 

We can predict the TV-dependence of D^ (TV) for a 
large range of TV values by combining the empirical 
TV-scaling [H HH (see inset Figg|) of D c . . m (N) with 
Dm(N). The Rouse scaling (D c . . m oc TV -1 ) is known to 
be violated and substantially higher exponents (> 1.5) 
have been reported [2(| HH HH . 

To determine Dm (TV) we have used Eqj5] and [H by ex- 
plicitly calculating the sum in Eq[4j Fig. 4 displays the 
predicted TDii(TV) along with the experimental data from 
Ref . [ll| . In agreement, both indicate a transition to a TV- 
independent regime beyond TV « 100. Further included 
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Figure 4: Self diffusivity of Li , Dl% = Dm + D c . . m with 
Dm from Eq[5] (□). The dotted line shows the prediction 
for a pure Rouse-type motion. Experimental [la] Dli for 
PEO:LiCF 3 S0 3 with EO:Li=20:l at T=342 K is included (a). 
The inset shows Dm and D c . . m individually. Dn values ob- 
tained from simulations are marked by crosses. 

in Fig[4] is an estimation of under the assumption 
that no entanglement effects arc present for D c o m by 
simply extending the scaling in the Rouse-regime to all 
N. As expected, the large-N behavior does not depend 
on the specific form of D c o m . However, one might spec- 
ulate that the appearance of the minimum in Dm{N) is 
indicative of entanglement effects. 

Interestingly, within the present approach the value 
of T3 can be estimated from experimental data. For 
very large N — Nl one has (see Eq.([Ta[>) r 3 = 
(N s b 2 /(6TrD Li \ N=NL )) 2 /(t 2 \ n =n s ) where b is the statis- 
tical segment length of the polymer, T2/TR ~ 2.2 and 
tr\n = n s can be estimated from D c o m « Dh\n=n s (at 



some small N = Ns) using the Rouse prediction. For the 
experimental data in FigfJ] this yields T3 ~ 80 ns. Via 
NMR experiments [2l[ local relaxation processes may be 
probed to extract information about t\. 

It is likely that the identified motional mechanisms are 
generally applicable to understand the ionic dynamics in 
different polymer electrolytes. Furthermore we checked 
that the presence of a small fraction of ions which are, 
temporarily, not bound to a polymer does not change the 
value of Dm by more than 10%. Of course, corrections 
will become relevant for larger ion concentration due to 
their mutual interaction. Note that our approach aver- 
ages over the different structural realizations and thus 
include, e.g., temporary complexation of an ion by two 
polymer chains. 

In summary, we have elucidated ion dynamics in poly- 
mer electrolytes by extracting microscopic properties 
from simulation and expressing them in analytical terms. 
This extends previous phenomenological approaches like 
the dynamic bond-percolation model, by assigning its key 
concept, i.e. the presence of renewal processes, a specific 
microscopic interpretation. For long chains the scaling 
relation Dn cx 1/^/ti2 • T3 is obtained which goes be- 
yond the expressions from the DBP model. In any event, 
for the transport of lithium ions fast transitions of ions 
between different chains are vital. Since the expression 
Dr,i(ri, T2,Tz) is now available one can estimate the pos- 
sible range of ionic mobilities for linear chain polymer 
electrolytes for all N. 
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